The Annals of Applied Probability 

2007, Vol. 17, No. 2, 421-439 

DOI: 10.1214/105051606000000781 

© Institute of Mathematical Statistics, 2007 
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Let Xi,X2,... be independent variables, each having a normal 
distribution with negative mean — j3 < and variance 1. We con- 
sider the partial sums S n = Xi + ••• + X n , with So ~ 0, and re- 
fer to the process {S„ : n > 0} as the Gaussian random walk. We 
present explicit expressions for the mean and variance of the maxi- 
mum M = max{S„ :n > 0}. These expressions are in terms of Tay- 
lor series about (3 = with coefficients that involve the Riemann 
zeta function. Our results extend Kingman's first-order approxima- 
tion [Proc. Symp. on Congestion Theory (1965) 137-169] of the mean 
for (3 I 0. We build upon the work of Chang and Peres [Ann. Probab. 
25 (1997) 787-802], and use Bateman's formulas on Lerch's transcen- 
dent and Euler-Maclaurin summation as key ingredients. 

1. Introduction. Let X\,X2, ■ ■ ■ be independent variables, each having 
a normal distribution with mean —j3 < and variance 1. We consider the 
partial sums S n = X\ + • ■ • + X n , with Sq = 0, and refer to the process 
{S n : n > 0} as the Gaussian random walk. In this paper we present explicit 
expressions for several characteristics of the distribution of the maximum 

(1.1) M = max{5 , TV :n>0}. 

The distribution of M plays an important role in several areas of applied 
probability. In queueing theory, it typically occurs in a regime called heavy 
traffic (see [2, 18, 20, 26]), in which the load is just below its critical level, 
and so the queue is only just stable with relatively large queue lengths 
and waiting times. For the limiting waiting time W = lim n ^ 00 W n , with 
Wi = and W n+ i = (W n + X n ) + (with x + := max{0,x}), it follows from 
Spitzer's random-walk identities that W is, in distribution, equal to M. In 
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the context of queues and heavy traffic, Kingman [20] was the first to observe 
the relevance of M in his 1965 paper. He noticed among other things: 

Despite the apparent simplicity of the problem, there does not seem to be 
an explicit expression even for EM..., but it is possible to give quite sharp 
inequalities and asymptotic results for small /3. 

Indeed, Kingman showed that, for (3 { 0, 

(1.2) EM = — — c + 0((3), c«0.58. 

2p 

Determining the tail distribution of M is tantamount to computing level 
crossing probabilities of the Gaussian random walk, that is, for x > 0, {M > 
x} = {t(cc) < oo}, where t(x) = inf{n > 1 : S n > x}. This level-crossings in- 
terpretation makes that the tail distribution of M is important in sequential 
analysis and risk theory. Chang and Peres [10] derived an exact expression 
(2.1) for the expected value of the first descending ladder height (actually, 
they consider the first ascending ladder height for the Gaussian random 
walk with positive drift), denoted as ESV_, with r_ = infjre > l:S n < 0}, 
which by the relation ES T _ = -/3/P(M = 0) (see [2], page 225) leads to an 
exact expression for P(M = 0). They present KS T _ as a Taylor series about 
(3 = with coefficients that involve the Riemann zeta function, a consid- 
erable achievement that generalizes first-order approximations of Siegmund 
[22] and second-order approximations of Chang [9]. 

Ladder heights fulfill an important role in probability theory, both in the 
exact analysis of random walks (see [2, 14]), and in the asymptotic analysis 
of boundary crossing problems [23]. In the latter quantity of interest 

is the limiting expected overshoot, defined as E(5y)/(2E£ r ), r = r(0), for 
(3 = 0. This quantity can be shown to be — £(l/2)/\/27r ~ 0.5826, with ((z) 
the Riemann zeta function. The same quantity arises in sequentially testing 
for the drift of a Brownian motion [11], corrected diffusion approximations 
[22], simulation of Brownian motion [3, 8], option pricing [6] and thermody- 
namics of a polymer chain [12]. These applications have in common that a 
Brownian motion is observed only at equidistant sampling points. As it turns 
out (Kingman [20] presents c as (2-n)- 1 / 2 E.^=i [v^v 7 ™ + V^ 17 ^) 2 ]' 1 ^, 
which by Euler-Maclaurin summation can be shown to be -(2tt)- 1 / 2 C(1/2). 
Similar relations are the topic of Problem 602 posed by Glasser and Boersma 
in [15]), the c in (1.2) is in fact — C(l/2)/\ // 27r, so Kingman, albeit in dis- 
guised form, related EM to the Riemann zeta function already in 1965. We 
shall extend Kingman's approximation (1.2) to an explicit expression for 
EM, in the same spirit as Chang and Peres extended the results of Sieg- 
mund [22] and Chang [9]. Moreover, we present a similar expression for the 
variance of M, to be denoted by VarM. The new expressions for EM and 
VarM both concern Taylor series about zero with coefficients that involve 
the Riemann zeta function. 
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The maximum and the first ladder height have been studied in the general 
setting of random walks with generally distributed increments (see [4, 9, 23]), 
the Gaussian random walk being a special case. For this general setting, Tay- 
lor series for the expected first ladder height and the expected maximum are 
presented in [4]. On a formal level, the results of Blanchet and Glynn [4] gen- 
eralize our results and those of Chang and Peres [10]. However, finding the 
coefficients of the formal description of the Taylor series in [4] is a nontrivial 
exercise and requires the expansion of a characteristic function and the nu- 
merical evaluation of an integral (see [4], Section 6, in which an outline for 
this numerical procedure is given). This does not lead to exact expressions 
for the coefficients as in [10] or as in the present paper. 

We first derive the Chang and Peres result [see (2.1) below] in our own 
fashion. Like Chang and Peres, we start from a Spitzer-type expression 
for P(M = 0), take its derivative with respect to (3, rewrite the deriva- 
tive in terms of the Riemann zeta function, and finally integrate to obtain 
(2.1). For rewriting the derivative, Chang and Peres built upon the 1905 
paper of Hardy [17] and present an analytic continuation of the function 
Li s (z) = S«S=i n~ s z n , known as the polylogarithm or Jonquieres function. 
They were probably unaware of the fact that Li s (z) is a special case of 
Lerch 's transcendent [see (2.4)], for which the matter of analytic continu- 
ation has been established in full generality by Bateman (and/or the staff 
of the Bateman Manuscript Project); see [13], Section 1.11(8) and (2.5). 
Hence, although Chang and Peres [10] give a separate proof, their Theo- 
rem 2.1 should be attributed to Bateman. 

Our derivation of (2.1) — that incorporates Bateman's formulas and an 
asymptotic determination of the integration constant — sets the stage for 
the derivation of the new explicit expressions for EM and VarM. As an 
aside, we obtain the following asymptotic results for (3 J. 0: 

(1.3) EM = ^ + C -^B + ^ + 0(p 2 ) 

and 

(1.4) 

where C(V 2 ) ~ -1-4604 and C(— 1/2) ~ -0.2079. In comparing (1.2) and 
(1.3), (1.3) contains an additional term \(3. This term, and — ^/? 2 in (1.4), 
follow from a rather intricate application of the Euler-Maclaurin summation 
formula. The error terms in both (1.3) and (1.4) will be replaced by Taylor 
series with coefficients that involve the Riemann zeta function. 
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1.1. Structure of the paper. We present our main results in the next sec- 
tion. Section 3 is devoted to an exposition of our derivation of the Chang 
and Peres result. The proofs of the new expressions for the mean and vari- 
ance of the maximum are given in Sections 4 and 5, respectively. The new 
expressions for the mean and variance of M are alternatives for their Spitzer- 
type counterparts. The latter tend to converge more slowly for a decreasing 
drift P, whereas the opposite holds for the new expressions. We investigate 
this difference in speed of convergence in Section 6. Concluding remarks are 
made in Section 7. 



2. Main results. We present three theorems. The first, on P(M = 0), 
is essentially due to Chang and Peres [10], but we give a separate proof in 
Section 3: 



Theorem 1 (Chang and Peres [10]). The probability that the maximum 
of the Gaussian random walk is zero satisfies 

(2 „ m -- 0) =v-^±^(=fJ}, 

/or < /? < 



Then, largely motivated by Chang and Peres, but taking our own ap- 
proach, we prove the next two theorems. 



Theorem 2. The expectation of the maximum of the Gaussian random 
walk satisfies 

(22)VM - 1 . c(i/2) i a 2 - c(-i/2-r) (-a 2 - 

(2.2) EM-- + w + -^ + -^X: r!(2r + l)(2r + 2) ^— 
for < p < 2^F. 

Theorem 3. The variance of the maximum of the Gaussian random 
walk satisfies 

4/3 2 4 H 24 

(2.3) 

2/3 3 ^ C(~3/2-r) (-^ 

V / 2^^^!(2r + l)(2r + 2)(2r + 3) V 2 



for < (3 < 20F. 
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The key ingredients for obtaining the above series are Euler-Maclaurin 
summation and a result on Lerch's transcendent. Lerch's transcendent is 
defined as the analytic continuation of the series 



(2.4) <f>{z,s,v) = J2(v + n)- s z n , 

n=0 

which converges for any real number f 7= 0, — 1, — 2,... if z and s are any 
complex numbers with either \z\ < 1, or \z\ = 1 and Re(s) > 1. Note that 
C(s) := <3?(l,s,l). We shall use the important result derived by Bateman 
[13], Section 1.11(8) [with C( s i v ) := ^(l) 5 ^) the Hurwitz zeta function]: 

(2.5) <&(z, a , v) = ^=-^(ln l/zf- 1 + z~ v f] C(s - r, v)^^, 

z v „ r! 

which holds for | lnz| < 2ir, s 7^ 1,2,3, ... , and v 7^ 0, — 1, —2, .... 

3. Proof of Theorem 1. Prom Spitzer's identity for random walks [24] 
we have 

(3.1) P(M = 0) = expj - f; ^P(5 n > 0) 1 = expj - £ Ip(-/V^) 1, 

I n=l J I n=l J 

with P(-) the standard normal distribution function 

(3.2) P(a) = -^= P e-W^dx. 



00 



The second equality in (3.1) follows from the normality of S n . 
With F defined by 

(3.3) F(/3) = £ ±-A= / e~U 2 > 2 dx, (3 > 0, 
we have 

_i 00 (l/2)/3 2 n _ p -(l/2)/3 2 00 (l/ 2 )/3 2 n 

F>((3) = ^=J2 7= = r= E , 

(3.4) 



2tt V 2 



Then by (2.5), when < \j? < 2vr, 



F>(3) _ -e-^ 2 \ T(l/2) (1 N-V 
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(3.5) 



+ e (l/2)/3 2 £ C 
r=0 



(-u/2)/? 2 r 



r\ 



-fV 



-1/2 



1 



EC, 



r=0 



1/2)^) 



2\r 



r! 



with denoting the Riemann zeta function. Thus, we get 



1 



1 



r=0 



1 



(-(l/2)/3 2 ) 



2\r 



< /3 < 



The series on the right-hand side of (3.6) converges uniformly in (3 S [0, (3q\ 
when < /?o < 2^/tt; see (6.3). Therefore, when we integrate the identity in 
(3.6) from to (3 < 2y/n, we may interchange the sum and integral at the 
right-hand side, and we get 



(3.7) 



r=0 



r!(2r + l) 



where L = lim m (F(/3) + ln(3). 

We shall show that L = — ^ In 2. To that end, we note that 



1 1 Z" 00 o 

F{0) = y±-L e~ u du 

±1 n Jy/(i/2)pn 



n=l 
00 



(3.8) 



1 ^ 1 / 2 
1 



|l n (l- e -(W) 



1 



^(l/2)^ 2 nVv^^vW^ 



ln/3 + -ln2 + o(l), 



as /? I 0. The function 



(3.9) 



g(y) ■-- 



1 / 2 



e~ u du-e- y 



y>o, 



decays exponentially as y — > 00, while flr(y) = 0(y~ 1/2 ), y j 0. It is then 
routine to show that 



(3.10) 



00 ;>oo 

f3 2 Y j g{\P 2 n)^\ g(y)dy, 

n=0 1/0 



1/32 
4 
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The latter integral can be evaluated as 

- 1 -= / e~ u du - e~ y ) dy 
y V Vvr Jy/y 



2 r°° 2 \ 
-= / e~ u du-e- y )lny 



roo / 2 1 

(3.11) -J Q [-j=---y- 1 / 2 -e~y + e-y)]nydy 



^ roo roo 

—= / y~ 1 ' 2 e~ y Iny dy — / e~ y lnydy 



/o •/ o 

-^r'(l/2) -r'(l) = -21n2, 



by Abramowitz and Stegun ([1], 6.3.1-4, 258). Hence, L = — ^ln2 indeed, 
and so it is shown that, for < (3 < 2y/n, we have 

( o 12) F(a) _ lnB l ln2 1 f, C(l/2-r)(-l/2r/^+i 
(3.12) F(/?)--ln/3--ln2-- 7 =^ -^-^ , 

which, by (3.1), completes the proof of Theorem 1. 

To recapitulate, we started from the Spitzer-type expression (3.1), rewrote 
its derivative (3.3) in terms of Lerch's transcendent (3.4), applied Bateman's 
formulas to obtain a Taylor series (3.5), integrated the Taylor series (3.7), 
and finally determined the integration constant L. 

Remark 4. The integration constant could have been determined from 
the relation P(M = 0) = —f3/ES T _ (with MS T _ the expected value of the 
first descending ladder height; see Section 1) and using the fact that ESV_ = 
— l/y/2 for /3 = 0, as proven by Spitzer [25]; see also [21]. Alternatively, 
one could use the first-order approximation in [18], that is, P(M = 0) = 
y/2p(\ + o(l)) as p i 0. The primary purpose of this section, however, is to 
set the stage for the next two sections, in which there is no other way of 
determining integration constants than to apply asymptotic methods. 

4. Proof of Theorem 2. From Spitzer's identity [24], we know that 

oo , oo , (l/2)/3 2 n x 

(4.1) EM = ^ -E(#) = £ - PP(-PV^) ■ 

n=i n n=A V2vrn / 

We then have 
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(4.2) 

Now we consider 
(4.3) 

We have 



1 | l_g C(l/2-r) 



°° 1 roo „ , 

G(I3) = Y^ e~ x / 2 dx. 



(4.4) 



n=i V 2vr 



-(l/2)/3 2 n 



"(l/2)/3 2 



z = e 



- (1/2) /3 2 



Then by (2.5), when ±/3 2 < 2vr, 



G"(/3) 
(4.5) 



-(l/2)/3 2 



2tt 



r(3/2) fi p2 



e -(i/2)p \ 2 

3/2 



-3/2 



+ e 



(1/2)/? 2 



Ec 

r=0 



(-(i/2)/? 2 r 



7" 



^ 2 



2\/2\ 
Therefore, we get 



r=0 



(-(V2)/3 2 ) r 



V 



(4.6) 



G'(/?) = -/r 3 -#(/?); 
ff(/3) = 



1 ^ C(-l/2-r)(-l/2)y _ 



r=0 



We note that -ff (/3) is well behaved in < (3 < and that 



(4.7) 



A. 

dp 



G(J3) 



2/3 2 



1 



■H(J3). 



By integration from to (3, we thus get 



(4.8) 



1 g, C(-V2-r)(-l/2)^ 2r+1 



r=0 



r!(2r + 1) 
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We shall show that 

(4.9) lim( G(£ ) -^)=4 

To that end, we use the Euler-Maclaurin summation formula (see [7], Sec- 
tion 3.6, pages 40-42) 

N r N 11 

£/(n) = / f(x)dx + -f(l) + -f(N) 

n=l 1 

m R 

( 4.10) +E i 7tTy(/ (2fc ~ 1) (^)-/ (2fc - 1) ( 1 )) 

- f N f (2m) (x) B2m{x ~ [x ^ dx 
h 1 [ ' (2m)! X ' 

where the B n {t) denote the Bernoulli polynomials, defined by 

(4 11) j£_ = ^B n {t)* 

e 2 -l ^ 



n=0 



and the 5 n = B n (0) denote the Bernoulli numbers. We apply (4.10) for 
m = 1, N — ► oo and 

1 Z" 00 2 1 

(4.12) fs(x) = -= e~ u du=:g(5x); 5 = -e 2 . 

V 71 " "'v 7 ^ 2 

Hence, 

N , ^ JV 



1 Z" 00 o , 

G N (e) = ^-= e-* 2 / 2 dx = J2fs( 
n=l v 2vr n=1 

(4.13) =j\(Sx)dx + ^g(S) + ^g(NS) 



+ l -B 2 {g'{N5)-g'{5))5- f" 5 2 g"(5x) Mx - [xJ) cix. 



2 -v, v , ^ v „ y i * v / 2 

Letting iV — > oo and noting that for g(y) = f^fy e ~ u2 ^ u there holds that 
5 1 <?" ~^ exponentially fast as y — > oo, we get 

/■OO 1 1 

G(e) = / g(5x) dx + - -B 2 g'{5)5 
(4.14) 1 
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Since |JB 2fc (x)| < B 2k for < x < 1 (see [1], 23.1.13, page 805), the last inte- 
gral at the right-hand side of (4.14) can be bounded by 

5 2 \g"(5x)\\B 2 dx = ^5 J \g"(y)\dy. 

We further get 

(4.16) g>(y) = --^ =1 /( y ) = ^J = (2 I/ + l)> . 

Therefore, we see that 

(4.17) 5g'(5) = 0O5 1 / 2 ), 5 f°° \g"(y)\ dy = -Sg'(S) = O^ 2 ). 
Furthermore, 

/oo roo roo / 1 /-oo \ 

g(5x)dx = 5~ 1 j s g(y)dy = 8- 1 j s {-= j e~ u du) dy 

(4.18) ^ ^ ^ s 

5- 1 !°°(-y= r e~ u ' du) dy - 5- 1 [ C4= f°° e~ u2 du)dy. 
vr J^/y J Jo VV 71 " Jy/y J 



Then from g(5) = \ + ©((J 1 / 2 ), we get 

(4.19) r 1 /* f 4= /°° e "" 2 du ) ^ = i + °( 5l/2 ) 

Jo V V 71 " J v 7 ^ ' ^ 



and 



(4.20) = -=y / e~ u du 



1 Z" 00 _ 2 \ 
— = / e u du) dy 



OO /-OO 1 1 

-/ y~y~ 1/2 -e-y dy 

Jo v 71 "^ 



1 Z" 00 1 

Therefore, 

/oo 11 
5 (fe) ( i 2; = --- + 0( ( 5 1 / 2 ), 5 10. 

It finally follows that 

(4.22) G(e) = ($-\ + 0{5 1 '- 2 )) + \ Q + 0(^ 2 )) ; 8 = \e\ 
and we obtain (4.9). It is thus concluded that 

U93) G(g) _ 1 1 1 f C(-l/2-r)(-l/2r 2r+1 

(4.23) G{(i)- w ----j=^ -^-^ /3 . 
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Combining (4.1), (4.2) and (4.23), we then obtain 

oo 



oo -(l/2)n/3 2 



1i \/2 



(4.24) 



7rn 

oo -(l/2)n/3 2 

Y- 

n=l ^ 
1 



n=l 



1 , l_g C(l/2-r)(-l/2)^ 2r 



r=0 







1_ 1 1 g C(-l/2-r)(-l/2)^ 2r+1 



2/3 2 4 



r=0 



r!(2r + 1) 



1 1 1 

2/3 + 4^ + y^F 



g C(l/2-r)(-l/2)^ 2r 



.r=0 



| - ((-l/2-r)(-l/2)^ 2r+2 ' 



r=0 



r!(2r + l) 



Splitting off the term with r = and replacing the summation index r 
1,2,... by r + 1, r = 0, 1, . . . in the first series in (4.24), we get 

(4.25) EM = L + «p + ^ + ^ £ CC-1/2 - r)(-l/2r/^ 



2/3 



2tt 



r=0 



r!(2r + l)(2r + 2) 



5. Proof of Theorem 3. From Spitzer's identity [24], we get 



(5.1) 



1, 



VarM=]T-E((S+) 2 ), 
, n 



n=l 



which, using the normality of S n , yields 

1 



VarM= 



(5.2) 



(xy/n- j3n) 2 e x /2 dx 



[ TI\/2tT Jfiy/n 

£ ( (/? 2 n + l)P(-/3v^) - -^=V^e-^/A , 



where the second equality in (5.2) follows from partial integration. We have 
established earlier [see (4.6)] that 

(5.3) ' £ = ' g C(-V2 - , ; )(-l/2) > 1 



n=l 



r=0 
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Therefore, it remains to evaluate 



(5.4) 



/(/?) = ]T nP(-/3v^) = E 



n 



71=1 



and to combine the results with (5.3) and (4.23) according to (5.2). 
There holds 



oo n 3/2 



-(l/2)/3 2 n 



-(l/2)/3 2 



(5.5) 



and by Bateman's result (2.5), 

. p -(l/2)/3 2 



-(l/2)/3 2 



3 

V = 1 



2tt 



r(5/2) ry^ 5/2 



"(l/2)/3 2 \2 



(5.6) 



+ e 



(1/2)/? 2 



Ec 



r=0 



(-(i/ 2 )/? 2 r 



r! 



3 (\f 



4v^V 
-3/T 5 



-5/2 x 



Ec 

r=0 



(-(1/2)/? 



2Nr 



1 - ((-3/2-r)(-l/2r f| 



r=0 



r! 



assuming that < (3 < 2^pii. The series on the last line of (5.6) is well be- 
haved in < < 2^^, whence l'((3) + 3/3~ 5 is integrable, and we obtain 



/(/?)- -r 4 = lim [1(e)- -re 



(5.7) 



10 V 



1_ g C(-3/2-r)(-l/2)^ 2r+1 



r=0 



r!(2r + 1) 



We shall show that 



(5.8) 



lim(/( £ )-fe- 4 ) = -^ 



by applying the Euler-Maclaurin summation formula (4.10) with m = 1, 
iV — > oo as before. We consider now 



(5.9) 



fs(x) 



Sx 



vr JJJx 



e u du=: h(5x); 
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in which 

1 f°° _ 2 

(5.10) h(x)=xg(x); g(x) = —= / e u du, x>0. 

V7T 



Then 



(5.11) 



1(e) 



h(Sx)dx + -/i(5) 



Next we shall take 5 J. 0, and to that end, we see that 



1 



1 1 r)V2 

(5.12) -h{5) = g(S) - -; h'(5) = g(S) - ^=e~ S ^ -, 

Furthermore, 

(5.13) - / h(5x)dx = -pr / h(x)dx — / h(x)dx, 
o Ji o z Jo o z Jo 

in which 

(5.14) 1 x 9 (x) dx i + 0(^/2). 

Also, by partial integration, 

/ h(x)dx= / / e _u (in I dx 

Jo Jo V^\J^ J 



610. 



(5.15) 



1 



e~ u du 



oo /"oo j.2 







o 2 



lx-V 2 e- x dx 



Therefore, 




(5.16) 


-f 

SJi 


Finally, 






h"(x) 



4^r Jo 

h(5x) dx 



oo 13 

x 3 / 2 e~ x dx = —^T(5/2) = 



40F 



16 



7 + 0(<5 1/2 )- 



165 2 4 



/i"(x) = (xg(x))" = 2g'(x) + xg"(x) 



(5.17) 

and 

(5.18) 



1 



x--)e- x eL 1 ([0,^)) 



1 



1 



2 5 2 (x- L xj)- 2 ^ 

k=l 



E 



cos 27rfcx 
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see [7], page 41. Therefore, 

poo 1 

6 h"{6x)-B 2 (x- [x\)dx 
Ji 2 

roc 1 

(5.19) =/ h"(x)-B 2 (x/S-\x/6\)dx 

Js 2 

1 1 



1 f°° 

^2~k?J s h "( x ) cos ( 2 ' K kx/5)dx^0, <5j0, 



27r2 , 
fc=i 

since J^ 00 h"(x) cos(2irkx / 5) dx — > as 5 [ by the Riemann-Lebesgue lemma 
on Fourier integrals. Putting this altogether, we find (recall <5 = ^e 2 ) 

(5.20) lim(/( £ ) - | £ - 4 ) = -1 + 1.1-1.1.1-0=-^. 
Hence, we obtain, for < (3 < 2y / 7r, 



(5 21) KB) 3 1 1 f- C(-3/2-r)(-l/2r ^. +1 

(5.21) 1W- 4 P 24 r!(2r + l) ^ ' 

We insert this result, together with (5.3) and (4.23), into (5.2) and get 

VarM = * t C(-3/2-0(-l/2r 2r+3 + 3 2 _ g 
V 7 ^^ r!(2r + l) ^ 4^ 24 

(5 - 2) H(27+l) * V "4 

1 ^C(-l/2-r)(-l/2)^ 2r+1 fl _ 2 



Splitting off the terms with r = and replacing the summation index r = 
1,2,... by r + 1, r = 0,l,... in the last two series in the right-hand side of 
(5.22), we get 

4/3 2 4 v^F 24 M 
1 » [ C (-3/2-r)(-l/2r 



^fo\ r!(2r + l) 

(r + l)!(2r + 3) 
(r + 1)! 



| C(-3/2-r)(-l/2r+^ 2r+3 



(5.23) + C(-3/2-r)(-l/r r 
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1 1 2C(-l/2) 1 2 

2 g((-3/2-r)(-l/2) r f+ 3 



2vr ^ r!(2r + l)(2r + 2)(2r + 3) ' 

6. Convergence comparison Spitzer formulas and Lerch series. It is im- 
mediately clear that the infinite series in (2.1), (2.2) and (2.3) converge more 
rapidly for smaller values of 0, while the contrary holds for their Spitzer-type 
counterparts (3.1), (4.1) and (5.1). To exemplify this difference in speed of 
convergence, we consider (4.2), that is, 

~ e -(i/2)n/? 1 » C(l/2-r) / 2 

The left-hand side series converges for all > 0, while the right-hand side se- 
ries converges for all EC, \0\ < 2^/tt. Prom [27] Section 13.151 (page 269), 

(6.2) 2 1 - s T{s)C{s) cos(isvr) = vr s C(l - s). 

With s = r + i, the asymptotics of the r-function and the fact that £(r + 
1/2) — ► 1 as r — > oo, we see that 

(0.3, ■4-C(l/* - -r)f_gy „ J__fg,, r^oo. 



2ir r\ V 27 vrV2r + 1 \47r / 

Hence, for comparing the convergence rates of the two series in (6.1), it is 
enough to find the point 0o > such that 

(6.4) e -d/2)/3 2 = 

47T 

With x = ^/3 2 , we need to solve xoe x ° = 2ir with xo > 0. This yields xq = 
1.4597, 0o = 1.7086, and the common value of the two members in (6.4) 
equals 0.2323. See [5], Section 2, where a similar strategy is developed in 
connection with the evaluation of Legendre's chi-function. 

7. Conclusions and outlook. We have presented analytic formulas of the 
Chang-Peres type, involving the Riemann zeta function, for the quantities 

OO -I 

(7.1) J fc (/?) = ]T-E((S+) fc ), & = 0,1,2, 

n=l 

yielding P(M = 0), EM and VarM of the maximum M of the standard 
Gaussian random walk with negative drift —0. The quantities can be ex- 
pressed using the normality of S n as 

<x> n k-l/2 f oo 
n=l 



(7.2) J k {0) = £ " / (y - 0f e -WW dy. 

±1 V2tt J (3 
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The following general result can be shown: There holds for k = 1,2, 

MP) 



(2/3)* 

k 



(7.3) + ± (*) ™^™ ( (^ _ \ 3 + ^.H-^ 

(-jjfc+ife! ~ ^(-fc-r + l/2)(-l/2) r ^ 2? - +fc+1 
+ \/2^ ^ r!(2r + l)---(2r + *; + l) ' 

when < /3 < 2y / 7r. That such a formula should hold for Jfc(/3) follows from 
differentiation of (7.2) A; + 1 times so that there results 



oo 



Jfc-1/2 



(7.4) Jf +1) (/?) = ("l) fc+1 fc! E ^e^ 2 ^ 



n=l 

in which the right-hand side is readily expressed in terms of Lerch's tran- 
scendent as we did before. Then, using (2.5) and integrating the identity 
thus obtained k + 1 times, we arrive at (7.3), except for the k + 1 integration 
constants that appear at the right-hand side of (7.3) as the coefficients of 
P 3 1 j = Oj 1) ■ • • >k. The actual determination of these integration constants 
is still based on the Euler-Maclaurin summation formula in its general form 
(4.10), but is so complicated that a full proof of (7.3) is outside the scope of 
the present paper. 

We are presently undertaking an effort to analyze a specific queueing 
model under a heavy traffic scaling. In the language of Queueing Theory we 
are dealing with a discrete-time bulk service queue with batch arrivals in 
which the arrival process is a Poisson process whose arrival rate A is just 
slightly smaller than the service capacity s. In the case that s = A + /3y/~\, 
with P > fixed and A — > oo (Halfin-Whitt regime, see [16, 18]), the equi- 
librium distribution of the queue converges to that of the Gaussian random 
walk. The analysis of this equilibrium distribution for finite A is, however, 
far more complicated than in the case of the Gaussian random walk. Already 
the first two moments require an evaluation in terms of Lerch's transcendent 
of the expressions 

(7.5) 5> fc_1/2 / e-W^xUx, 

71=1 13 

for all integer k and all z = 0, 1, The results obtained by us for the expres- 
sions (7.5) allow us to fully establish (7.3), that is, including the integration 
constants. 

A further comment concerns the restricted validity of formulas (2.1)- 
(2.3), namely for < f3 < while the corresponding Spitzer formulas 
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(3.1), (4.1) and (5.1) make sense for all /3 > 0. By using (6.2), the explicit 
formula £(s) = Y^n=i n ~ s valid for s > 1, and some further manipulations, 
the infinite series occurring in (2.1)~(2.3) can be re-expressed as follows. We 
have for (3 > 

C(l/2-r)(-(l/2)/? 2 r 



(7.6) 



(7.7) 







E 

r=0 



r!(2r + 1) 



v/2¥ vr 



*/ 4 S 



/* 2 



E 

r=0 



2vr 2 



Re 



i/4 5i 



4-7T 



-if 
47T 



2\r 



r!(2r + l)(2r + 2) 



and 



(7.8) 



in which 



(7.9) 



(7.10) 



(7.11) 



2(3 3 



E 



2\r 



C(-3/2-r)(-(l/2)/? 
2vr r!(2r + l)(2r + 2)(2r + 3) 



P 3 

47T 3 



Re 



4-7T 



So(6) 



Si (6) 



5 2 (6) 



(arcsin(fe/n) 1 ^ 2 



(6/n 



n=l 



'TT 

2~) 



E 

n=l 



n 



n 



Vn-b), 



'TT 

46 



E ~^(Vn- Vn-b) 



n=l 



are well defined functions for 6 (= — i(3 2 /Att) in an open set containing the 
imaginary axis. (These alternative expressions are intimately related with 
Lerch's transformation formula ([13], 1.11(7) on page 29) and complement 
the results obtained through analytic continuation by Chang and Peres [10].) 
Thus, for values of (3 larger than 2^/tt, we can evaluate P(M = 0), EM, 
Var M using (2.1)-(2.3) and (7.6)-(7.11). Although the series in (7.9)^(7.11) 
converge slowly, they can be evaluated quite conveniently by using dedicated 
forms of the Euler-Maclaurin summation formula (4.10). We may also note 
that the infinite series occurring in the expression (7.3) for the general Jk(/3) 
can be re-expressed in a similar fashion. 
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Finally, let us return once more to the rather special constant —£(1/2) /v2~7r ~ 
0.5826 (see Section 1). From Theorem 2 [or even (1.2)], we get 



where 1/(2/3) is known to be an upper bound on EM that becomes tight in 
heavy traffic, that is, for (3 J, (see [19, 20]). Hence, the difference between the 
upper bound and the true value tends to a nontrivial constant —£(1/2) / \[2tx . 
In [12], constants of this type are studied for a general class of random walks 
(possibly non-Gaussian) with zero drift. From Theorem 3 and the general 
result (7.3), we see that such nontrivial constants exist for higher moments 
of M as well. 
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